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ABSTRACT 

Analysis and Evaluation of a Forward-Viewing Scanning 
Radar Scatterome ter System. (August 1970) 

Richard Harrison Arnold, B S., 

University of Southwestern Louisiana 
Directed by Dr. J. W. Rouse, Jr 

A scatterome ter is a radar system used to measure 
the radar cross-section of area extensive targets, par- 
ticularly terrain As a part of the evolvement of this 
instrument, it has been proposed by a leading aerospace 
company that a scanning scatterome ter be developed which 
would integrate for 0 2 seconds the returning signal from 
a terrain element defined by a minimum doppler frequency 
difference of 280 Hz. This system would provide measure- 
ments of the scattering coefficient of an area ±35 de- 
grees transverse to the aircraft direction The purpose 
of this research effort was (1) to determine if the al- 
lowed integration time was sufficient for a given terrain 
element size for an accurate determination of the scat- 
tering coefficient, and if necessary, (2) to develop 
methods for improving the measurement accuracy without 
altering the basic scanning technique. 

A scanning s cat ter ome ter is analyzed by system 
subassemblies transmission, reception, signal reduction 
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and logic, recording, and CRT display. Possible klystron 
transmission frequency drift problems and potential cali- 
bration problems are discussed. The calibration system 
is analyzed with respect to system operations An al- 
ternative system is presented and shown to be equal to or 
better than the proposed calibration system. 

The accuracy of the scattering coefficient deter- 
mination was shown to be directly related to the signal 
integration time and the terrain element size. The pro- 
posed integration time for the terrain element specified 
by the given doppler frequency difference was proven to 
be insufficient. Either an increase in the system's 
integration time or the element size is proven necessary 
to have sufficient accuracy. 

An analysis of the possible use of frequency agil- 
ity by the scatterometer system showed that definite ad- 
vantages do exist. The use of frequency agility is 
shown to increase the accuracy of the scatterometer 
without increasing the integration time or the terrain 
element size. Further analysis involving defining a more 
exact analysis of the effect of frequency agility upon 
the accuracy of the scatterometer is suggested 
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CHAPTER I 
INTRODUCTION 

A scatterometer is a radar system used to measure 
the radar cross-section of area extensive targets, par- 
ticularly terrain. As a part of the evolvement of this 
instrument, it has been proposed that a scanning scat- 
terometer be developed [1] This system would provide 
measurements of the scattering coefficient of an area 
±35 degrees transverse to the aircraft direction How- 
ever, the engineering requirements for accomplishing the 
scanning operation demands a strict limitation on the 
time allowed for receiving the backs cattered energy from 
each terrain element The purpose of this research ef- 
fort was (1) to determine if the allowed integration 
time was sufficient, given a specific terrain element 
size, for an accurate determination of the scattering 
coefficient and, if necessary, (2) to develop methods 
for improving the measurement accuracy without altering 
the basic scanning technique 

The citations on the folloitfing pages follow the 
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The Scatterometer 

Radar systems belong to the class of instruments 
termed remote sensors. These instruments are utilized 
to detect characteristics of objects from a distance m 
somewhat the same manner that the eye observes objects 
at a distance Remote sensors have been found to be use- 
ful, especially in earth resources studies, since they 
have the capability to observe objects under conditions 
where physical contact is impossible [2] 

Radar systems are termed active sensors because they 
provide their own source of illumination This illumina- 
tion is accomplished by transmitting electromagnetic 
waves Due to the nature of these waves , reflection 
from objects occurs and a portion of the reflected trans- 
mitted signal is directed towards the transmitter 

Radar systems are capable of determining factors 
such as time delay, frequency difference, and power loss 
between transmission and reception From these three 
factors, certain characteristics of targets can be de- 
termined. Time difference can be used to determine the 
distance between the radar transmission source and the 
reflecting object. Frequency difference can be used for 
deteiminmg relative velocity between the source and the 
object, or it can be used for spatial location using the 
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principle of the doppler frequency shift 

Radar return is the result of signal reflection from 
an object Most objects, especially terrain surfaces, 
which radar systems are utilized to detect, have what is 
considered as "rough" surfaces Upon striking a rough 
surface, the radar signal will be dispersed m some 
pattern dependent upon the surface configuration For a 
large surface, such as the ground, many points within the 
surface constitute the object, or scatterer. For these 
large surfaces, radar return is subject to the scattering 
pattern of the sum of all the individual scatterers 
With a slight change m relative position between the 
radar source and the object composed of many individual 
scatterers, the scattering pattern may change consider- 
ably This causes the signal amplitude of the radar re- 
turn to become random when motion occurs between the 
source and the object [3] 

Parameters, other than the scattering pattern, af- 
fect radar return Some of these are angle of incidence, 
polarization, and transmission wavelength. These para- 
meters are functions of the radar transmission Surface 
parameters affecting radar return are, m addition to 
roughness, electrical conductivity and permittivity All 
of these parameters will affect the power level of the 
radar return Those parameters which are functions of 
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radar transmission can be varied to observe power varia- 
tion as a function of the parameter variation. Some- 
times, these observations are useful m determining prop- 
erties of the reflecting object One of these, variation 
of the angle of incidence versus a type of normalized 
power loss called the backs catter coefficient, is the 
mam objective of the specific mode of radar system 
called the s catterometer . 

A radar m the s catterometer mode measures the re- 
turn power received from a ground area or object which 
scattered the signal transmitted by the radar. As 
stated, the primary objective of this radar mode is to 
obtain values of the radar backscatter coefficient, also 
called the differential scattering cross section, for a 
particular surface at a particular angle of incidence. 

As with other types of radar systems, the s catterometer 
can be designed to transmit and/or receive with different 
polarization schemes and to transmit at any particular 
microwave wavelength. Scatterome ters are generally de- 
signed to determine the backscatter coefficient at sev- 
eral angles of incidence An ordinary radar system is 
capable of determining the backscatter coefficient for 
many types of objects but generally at only the parti- 
cular angle of incidence at which the object is illumi- 
nated, and determining the angular, wavelength, or polar- 
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ization dependence of the backscatter coefficient is not 
its primary objective. A s catterome ter has become a 
distinct radar mode due to its peculiar type of opera- 
tion 

A surface is either rough or smooth relative to the 
incident wavelength A particular wavelength will de- 
termine the degree of roughness of the surface. This de- 
gree of roughness has the dominant effect upon the mag- 
nitude of the backscatter coefficient Since the degree 
of roughness will determine to a certain degree the 
scattered field pattern for terrain, the scattered field 
energy measured at any incidence or aspect angle is a 
function of the terrain roughness The large-scale ter- 
rain features will determine the magnitude of the back- 
scatter coefficient for a small incidence angle, and for 
a large incidence angle, the small-scale structure de- 
termines the magnitude [4] 

The overall effect of either the polarization of the 
incidence wave or the complex dielectric constant has 
not been adaquately determined For studies of specific 
materials, where their effect has been empirically de- 
termined, their utilization can prove effective [5] 

Scatterome ter data has proven to be useful m de- 
termining differences m sea ice types and sea state con- 
ditions The scatterome ter has been shown capable of 
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determining the roughness variations m the homogeneous 
surfaces of the sea and ice. Scatterometer data taken 
from surfaces having non-homogeneous materials has been 
largely unsuccessful m determining the surface charac- 
teristics 

Scatterometer data taken over the Arctic ice cap has 
provided backs catter coefficient versus incidence angle 
plots which have been used to distinguish young ice from 
old ice [6 ] . For the past several years, scatterometer 
data concerning sea reflectance has been collected. The 
assumption made to justify scatterometer utilization is 
that as the seas become more turbulent, the water sur- 
face becomes rougher. Investigations have shown that a 
correlation exists between sea state conditions and the 
backscatter coefficient [7] . However, recent work has 
also shown that there apparently exists a sea turbulence 
value above which the backscatter coefficient becomes 
useless for determining exact sea states [8] 

Numerous investigations employing radar m the 
scatterometer mode are presently in progress. Texas A§M 
University [9], University of Kansas [10], Naval Research 
Laboratories [11], Ohio State University [12], Waterways 
Experiment Station [13] , among others are attempting to 
further refine the capabilities and applications of this 


sensor. 
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The Scanning Scatterometer 

Many laboratory experiments have been made m order 
to obtain backs catter coefficient versus incidence angle 
plots for various soils, road types, and crop types 
Some of this data has been collected to provide a catalog 
standard for scatterometer data analysis [14] 

Present airborne scatterometer systems view beneath 
the aircraft within a very narrow beamwidth For exam- 
ple, the NASA MSC scatterometer views beneath the air- 
craft within a ±1 5 degree beamwidth This drastically 
limits the terrain coverage per flight. Therefore, a 
more efficient viewing scatterometer, one which is capa- 
ble of viewing more surface area per flight, is desir- 
able Such a scatterometer is the proposed scanning 
scatterometer 

The scanning scatterometer is designed to view 
ground area transverse to the aircraft as well as along 
the flight direction. This is accomplished partially 
by sweeping the scatterometer antenna 1 pattern coverage 
transverse to the aircraft direction Basic to the 
scanning scatterometer is the use of the doppler fre- 
quency shift to help delineate the specific ground loca- 
tion to be viewed No known sufficient quantitative 
analysis of the ability of the scanning scatterometer to 
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effectively determine backs catter coefficient values has 
been made. This effectiveness is directly related to the 
use of the doppler frequency shift. Therefore, such an 
analysis will be partially concerned with the effect of 
using the doppler frequency shift This research de- 
velops this analysis for a specific scanning scattero- 
meter to be described m Chapter II. 

Summary 

In this document, a scanning scatterometer system 
is described using block diagram type notation. The 
system is presented m sections associated with the scat- 
terometer subassemblies, transmission, reception, signal 
reduction and logic, recording, and CRT display, in 
Chapter II. This Chapter is the basis for subsequent 
examinations m Chapters III and IV. 

A theoretical analysis of the scanning scatterometer 
system was performed m Chapter III with respect to the 
determmance of the backscatter coefficient. Due to the 
stochastic nature of a backscattered signal, this anal- 
ysis was conducted using probability theory. The signal 
power spectrum is traced through the signal reduction 
section The resultant signal spectrum after half-wave 
linear detection is computed. The number of independent 
samples available during the integration period was de- 
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termmed From this information, the accuracy of the 
ability of the scatterometer to determine the backs catter 
coefficient was calculated by using the average power 
variance and the average voltage squared variance 

The advantages of implementing frequency agility on 
future scanning scatterometer systems are determined m 
Chapter IV with respect to previous studies [15] con- 
cerning frequency jumping radar systems The loss of 
accuracy m determining the backs catter coefficient by 
frequency agility was weighed against the possibility of 
increasing the number of independent samples available 
during the integration time The results show that fre- 
quency agility implementation is advantageous from the 
standpoint of increasing the accuracy of the backscatter 
coefficient determinations 

A summary of the study and recommendations for 

further work is presented m Chapter V. Further work is 
< 

recommended to determine the full potential of using 
frequency agility on future scatterometer systems. 
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CHAPTER II 
SYSTEM 

Transmission. 

Transmission, as with any radar system, is described 
m terms of transmitting and receiving periods, signal 
features, and antenna pattern This system is designed 
with the intent to incorporate two altitude ranges, thus 
requiring dual transmission modes. Therefore, transmis- 
sion and power specifications will be discussed separate- 
ly from the antenna design and pattern. 

The transmitted frequency is the interrupted con- 
tinuous wave (ICW) type The system operates at 13 3 GHz 
using a klystron power output of ten watts minimum The 
pulse repetition frequency (PRF) of this signal is a 
function of aircraft velocity and altitude m each alti- 
tude range. This is also the case for the transmitting 
and receiving periods (see Fig. II-l) [16] 

Since ground return is desired from points up to 
50 degrees from the aircraft nadir, reception time must 
at least equal twice the time for a signal to be trans- 
mitted to the furthest point Therefore, 

*T* — ^ farnox 

! r C COS-0 


( 2 - 1 ) 
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Knowing the transmitting and receiving periods, the 
pulse repetition frequency is found to be 

+ 2 T d 



Substituting the value for maximum receiving time 

i _ + 2 ~Q) GO 5 O h ^ hrrax 

C COS 0 (2-2) 


where C = speed of light 
0 = incidence angle 
~"Q = transmitting time 

T, = time delay between transmission and re- 
d ception 

maximum aircraft height 

This PRF is subjected to a lower limit due to the 
system. The transmitted signal is subjected to a fre- 
quency shift upon striking the ground This doppler fre- 
quency shift (£^) is used in the system to delineate 
ground area, hence by the sampling theorem, the minimum 
PRF is 


PRF mn 



mat 


(2-3) 
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For a view at a constant angle of incidence, a maxi- 
mum f^ occurs at the ground area m the flight direction. 
Doppler frequency shift along this line is defined as 

f , - sin 9 = sin 0 

A c (2-4) 

Therefore, taking theta maximum as 50 degrees, and 
maximum velocity as 300 knots (154 meters per second) 

f = E(l ^ Xa3 --- O ' -5in0 = 10.17 KHz 

« mot* 3X| 0 * (2 _ 5) 

311,1 SO. IhKHz (2-6) 

where 7/= aircraft velocity 

"F = transmitted frequency 
X- transmitted wavelength 
C = speed of light 

Two altitude ranges are selected of which the mini- 
mums were 2,000 ft. and 5,000 ft Using a time of one 
microsecond for T^, two transmission modes are obtained 
and shown m Table II-l. 

The signal is continuous when emitted from the kly- 
stron. Pulse modulation is accomplished by the use of a 
duplexing switch, the maximum switching rate being 45 KHz. 



Altitude 

Transmitting 

Pulse 

Pulse Rate 
Frequency 

2,000 - 
5,000 ft. 

4.08 usee 

45 KHz 

4,000 - 
10,000 ft. 

8.16 ysec. 

22.5 KHz 


TABLE II-l TRANSMISSION MODES 


•p* 
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The scatterometer uses a klystron m the trans- 
mission section (see Fig. II-2) The klystron is a high 
frequency signal generator and amplifier. The reflex 
klystron or a multicavity klystron with feedback is most 
commonly used for this purpose Klystrons, being com- 
pact, efficient at high frequencies, and of simple con- 
struction, are utilized m place of tubes. Tubes gener- 
ally have relatively large electrodes and long lead wires 
which create capacitance and inductances too 'large for a 
microwave resonant circuit Also, the electron transit 
angle is large for tubes at microwave frequencies causing 
an output delay m the tube (i e. increased transadmit- 
tance) . To minimize these effects, tubes can be mini- 
turized, but this limits power-handling capabilities [17]. 

The klystron circumvents these problems. First of 
all, it can be of compact construction. Second, it uti- 
lizes cavities for the microwave resonant frequency, and 
third, the electron transit angle can be minimized 

The klystron (Fig. II-3) is basically a tube with 
cavities operating on the principle of electron velocity 
modulation. Amplification occurs by feeding a signal 
into the input cavity, called the buncher At any one 
time, the voltage across the buncher gap will be v = V 
sm yt, where V is the input voltage and y is the signal 
frequency Electrons passing the gap when v^ = V will 
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travel slower than electrons passing the gap when v = 0 
Electrons will travel their fastest when v = V Thus, 
the electrons are velocity modulated by the voltage 
across , the buncher gap 

These electrons are emitted from a cathode and fo- 
cused into a beam by the electron gun. The cavities are 
at a positive potential for beam acceleration and not for 
electron attraction purposes After the electrons are 
velocity modulated by the initial cavity, they tend to 
bunch up m the drift space Electrons accelerated dur- 
ing the high voltage period overtake those slowed by the 
previous low voltage period, thus causing the bunching 
effect. They then pass through the catcher cavity and 
collide with the collector 

The catcher cavity, or output cavity is placed such 
that its gap is at the place of maximum bunching Here, 
power is extracted from the bunched electrons The en- 
ergy received from the bunched electrons is generally 
many times greater than that required to veloci ty-modu- 
late the electron beam. Thus, amplification occurs 
Oscillation can be obtained very easily by coupling the 
buncher and catcher with a feedback circuit. The fre- 
quency will be the resonance frequency of the oscillator 
having a no-load condition 
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In multi cavity klystrons, several buncher cavities 
are used. Proper adjustment of the phase and magnitude 
of the waves m these bunchers produces a very tightly 
bunched electron beam. This additional bunching in- 
creases the output gam. 

A reflex klystron (Fig. II- 4) utilizes only one 
cavity to achieve oscillation. Electrons are emitted by 
the cathode and are velocity modulated when initially 
passing through the cavity grids These grids are held 
at a positive d-c potential with respect to the cathode 
and will have the voltage operation shown m Fig. II-5. 

By proper adjustment of the beam acceleration and 
repeller voltages, the variation of voltage on the grids, 
caused by initial oscillation, will cause resonance. The 
electrons are repelled by the highly negative reflector 
and pass back through the grid gap. If the bunched elec- 
trons are slowed down by the gap-voltage cycle, their en- 
ergy is applied to the oscillation. 

The duplexer switches are devices which switch an 
antenna either to the transmitter or receiver. These 
switches are pieces of ferrite strips whose impedance is 
controlled by an externally applied d-c magnetic field 
Signals create magnetic fields which, using the non- 
isotropic property of ferrite, create low loss trans- 
mission m one direction and almost total power absorp- 
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tion m the other direction. The switches are designed 
such that the maximum switching rate is 45 KHz. When one 
switch is utilized to terminate flow of power from the 
klystron to the antennas, another switch allows passage 
of the signal reception to the receiving section and 
vice-versa Thus, the klystron is allowed continuous 
oscillation while the scatterometer operates as an ICW 
system. 

The scatterometer output power can be transmitted 
either horizontally, vertically or time-shared between 
the two types to afford selective surface scattering 
modes This is accomplished using another ferrite switch 
which accomplishes the task of a ferrite circulator. 

Two antennas are utilized to provide adequate iso- 
lation and dual polarization, one horizontally polarized 
and the other vertically polarized The antenna is a 
slotted waveguide m a duplexed scanning phased array. 

The beamwidth is about 120 degrees fore-aft and 3 degrees 
port-starboard Each hole on the slotted waveguide is 
associated with a ferrite phase shifter A beam steering 
computer is utilized to vary the current m the ferrite 
to change the phase such that the desired lobe position 
is obtained. The computer is designed to allow for 27 
angles m the transverse direction for the lobe. The 
rate of transverse angular change is controlled by the 
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factor v/h, where h is the height of the aircraft and v 
is the velocity of the aircraft. 

The ground area covered is dependent upon the time 
length spent viewing at each angle Certainly, total 
ground coverage is impossible using this scanning proce- 
dure, however, this is not the purpose of this design 
Instead, it is a general idea of ground roughness type, 
m the form of the backscatter coefficient, for compar- 
ison with coincident microwave radiometer data This 
matter will therefore have little significance m the 
equipment design. The ground coverage will be as shown 
m Fig. II-6. 

Within this ground coverage, total radar backscatter 
will be placed on tape, however, the visual display will 
utilize only that portion of the data corresponding to 
the radiometer area coverage To do this, only data re- 
turning from a specified ground area must be filtered 

/ 

out of the total data return The method of filtering 
will be discussed later The ground area defining this 
specific area is referred to as the cell The cell is 
defined as having 75 foot length m the direction of 
flight (thus defining the doppler bandwidth) and a 3 de- 
gree areal coverage across Also specified for thi,s 
defining cell is the aircraft height of 1,000 foot alti- 
tude and viewing angle of 50 degrees from the vertical m 
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FIG I I- 5 REFLEX KLYSTRON VOLTAGE PATTERN 
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FIG I I- 6 SCATTEROMETER GROUND COVERAGE 
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the flight direction ((f) = 0). This determines the center 
frequency and doppler bandwidth of the area. As the 
fan beam angle moves off of (j> = 0 , a variable low pass 
filter is adjusted such that the area viewed along a 
circular path 50 degrees from the vertical has the same 
ground length. The angle of scanning for the CRT data 
is always 50 degrees from the vertical, and will make a 
circular scan of the ground sweeping a ground arc of 72 
degrees relative to the aircraft (Fig. II- 7) 

The cell definition was designed to make the scat- 
terometer and radiometer viewing screens complementary. 
Even so, the basic cell can be adjusted if a more accu- 
rate determinance of backs catter coefficient is found to 
be necessary To do so, the fact that all ground cover- 
age depends on the definition of the basic cell must be 
considered to see if system limitations are exceeded. 

Re cep t ion 

The returning signal is a function of the terrain 
scattering coefficient. It has a center frequency at 
the transmitted frequency 13.3 GHz but contains an en- 
velope of the doppler frequency shift (f^) ♦ In measuring 
the backscatter coefficient, f^ is utilized to specify 
ground position via frequency The envelope of the f^ 
spread is somewhat as shown m Fig II -8 The backscatter 
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from vertical incidence is much greater than that from 
angles away from the vertical. When the fan-beam antenna 
lobe is pointed off vertical, the center frequency spike 
tends to demmish m size. 

The frequency spread (approximately 21 KHz) will be 
constant m relation to the center frequency regardless 
of mixing. The spread in percentage of the center fre- 
quency, will change. 

4 spread - X3. 3 GHz ” 1 ’ 6 x 10 ' n 

% spread - J* -g § f “ 8.4% 

Increasing the percent of spread about the center fre- 
quency allows for less stringent design specifications 
for signal filters, such as the trap and low-pass filters 
Another consideration concerning the return signal 
is that of conserving the return signal power Evalua- 
tion of the backscatter coefficient (oo) is directly de- 
pendent upon the magnitude of return power from a ground 
area. Since the determmance of o 0 will have a variance, 
any uncorrected loss of power will have an adverse ef- 
fect upon the variance. Therefore, steps must be taken 
to conserve power. Taking these two facts into con- ' 
sideration, the returning signal is subjected to fre- 
quency mixing to reduce the center frequency and is sub- 
jected to a calibration system to eliminate amplification 
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fluxuations . This affords easier signal handling and de- 
tects hardware effects upon power. The reception section 
of the s catterometer system is shown m Fig. II-9 

The signal is received by the polarized antennas and 

t 

passes through a circulator to the reception switch. 

When signal reception is desired, the signal proceeds 
through the mixer, where the center frequency is reduced, 
and to the logarithmic amplifier. The logarithmic ampli- 
fier is designed to equally amplify all portions, or fre- 
quency amplitudes. The trap filter is a band-e limmator 
filter which eliminates the high amplitude center fre- 
quency of the return signal. This eliminates the possi- 
bility that other circuit amplifiers will saturate. The 
problems with the mixer circuit are of great interest 
and will be discussed below. 

The principle of frequency changing is basically 
simple. To accomplish the frequency change, the re- 
turning signal must be mixed with a signal at a frequency 
equal to the sum of the original center frequency plus 
the desired center frequency. The type mixer used and 
its construction is shown m Fig. II -10. 

The mixing frequency is constructed such that the 
center frequency of the return signal will be added to 
the local oscillation of 250 KHz. To accomplish this, 
certain assumptions must be made. The klystron oscilla- 




SSBM = Single Sideband V.A. = Variable Attenuator 

Modulator 


FIG II -9 SCATTEROMETER RECEPTION SECTION 


to 
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FIG 11-10 SINGLE-ENDED MIXER 



13.3 GHz 


Sf = variation about the center frequency 
FIG. II- 11 EFFECT OF FREQUENCY TRANSMISSION DRIFT 
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tion at 13 3 GHz is subject to drift over a small range 
due to cavity geometry This will appear as shown m 
Fig. 11-11 This drift characteristic will vary with 
emitter and catcher voltage changes The power monitor 
provides constant amplitude, yet the drift is certain to 
occur If a drift m frequency has occurred since the 
transmission and mixer reception of the signal, the fre- 
quency spectrum after mixing has occurred will be 


(fc * & f c + 250 Kttz ) ~ - £ foj 

= 2 50 KHz + ££ ± (2-7 ) 

where <C = variation 

= center frequency 
= doppler frequency 

Any <5f c will cause the signal reduction circuit, 
using doppler frequencies, to locate an area other than 
the desired area The possibilities of this happening 
can be avoided. 

The maximum travel time, t , for the transmitted 
signal is 


2h may . 2(10") 1 

' c cos o ' costo” “fya-fo* 
— . 3 2 X I 0* Seconds 


(2-8) 
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where C = speed of light 

h = aircraft altitude 

The doppler frequency change (Af^) m a cell block, 
say Af^ - 287 Hz, is chosen to allow for an adequate 
number of independent samples to be received Therefore, 
for the cell to be completely missed m analyzation 


Afj = £f e * 287 Hi 


The frequency change per second oc curing during the 
transmission time will be 


K = 

i~ 

‘-'max 


2 £7 

.32x/o'‘ ) 


IOHUz/sec 


(2-10) 


Klystrons having a frequency time drift of one-one - 
hundreth of this magnitude are available Therefore, 
proper selection of the ground area should not constitute 
a problem. 

The local oscillation is introduced with the use of 
a single sideband modulator (SSBM) This SSBM receives 
the signals from the klystron and a local oscillator 
The local oscillation frequency of 250 KHz modulates the 
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13.3 GHz klystron signal creating signals of center fre- 
quency with upper and lower sidebands 

The carrier frequency is denoted by 

X(t)= Pi Cos GJ k t (2-11) 

where X(t) = carrier signal 

is klystron frequency 
Q — is amplitude 

The amplitude. A, of the signal after single sideband 
modulation is 


P c + k, CO 5 U> m t 


( 2 - 12 ] 


where. 03 ^= is modulating frequency 

/? c = is a constant value 

k, = is a constant 

The signal after single sideband modulation now becomes 

y(t) 

kjflm cos cos oO k t 


-I- CO$(GO k +-oO„)-t 


A 




(2-13) 
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where fp a 


K A- 

a 


The resultant frequency spectrum will appear as 
shown m Fig 11-12 

If there exists more than one modulating frequency, 
the spectrum will appear as shown m Fig 11-13 

Through all the reception hardware, the return power 
must be conserved and its amplication specifically known 
‘to eliminate any increase m Oo variance To accomplish 
this, a process of signal calibration is used to notice 
and correct for any system effect on power The system 
is designed to allow for a 4 0 db maximum loss m the 
system and for a 8.0 db maximum system noise level at 
250 KHz If this is allowed to affect the variance, it 
will destroy the effectiveness of the system The ef- 
fects of a slight power loss can destroy the system's 
ability to determine crop types, as can be seen from the 

crop db differences m Table II-2 [14]. 

\ 

Table II-2 data implies that a system having a two 
decibel variance m the o 0 calculation can not determine 
between several crop types. The ability to construct 
the entire c 0 vs 0 curves is present, thus allowing 
identification by curve signature instead of only one 
angle. The problem of justifying the use of a single an- 
gle for identification depends upon eliminating the vari- 
ance subjected to o 0 by the system. 
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. 11-12 MODULATION FREQUENCY SPECTRUM 



FIG 11-13 

MULTIPLE MODULATION FREQUENCY SPECTRUM 
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TABLE II -2 

BACKS CATTER COEFFICIENT DATA AT A 50 DEGREE 
INCIDENCE ANGLE AND VARIOUS TRANSMISSION FREQUENCIES 


a 0 Power 


Material 

Level In 
Negative db* 

Frequency 

Used 

Weeds on Dry Terrain 

18 

db 

34 5 

GHz 

Weeds Over Wet Terrain 

25 

db 

9 4 

GHz 

Wet Soybean Stubble 4 Feet Tall 

10 

db 

35 

GHz 

Green Soybeans 3 Feet High 

14 

db 

35 

GHz 

Oats 1/2 Inch High 

16 

db 

35 

GHz 

Oats 4 Feet High 

14 

db 

15 5 

GHz 

Wheat Stubble 12 Feet High 

12 

db 

35 

GHz 

Green Wheat 1/2 Inch High 

16 

db 

35 

GHz 


The transmitted power was used as the reference for the 
db scale 
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The method used to eliminate this variance is to 
inject another signal into the system and allow it to be 
subjected to all the return signals manipulations. The 
final calibration signal is then divided into the end 
results of the manipulated return signal to remove all 
gain and noise variations . 

One proposed method of inserting this calibration 
signal is shown by the hardware schematic Here the 
calibrated signal combines with the local oscillator 
signal prior to entering the SSBM The SSBM amplitude 
modulates the klystron generated signal with the input 
signals . 

The carrier and lower sidebands are eliminated and 
only the upper two sidebands emerge from the SSBM. This 
signal is injected into the mixer to produce a reduction 
m frequency of the returning signal. Both signals are 
equally processed. The final calculation of power for 
the 265 KHz signal is divided into the terrain signal 
power This corrects for gam and noise variations from 
the mixer to the divider. 

Several problems are encountered by utilizing this 
method The major difficulty is that the returning 
power will adjust about 250 KHz. This can most easily 
be shown graphically (see Fig 11-14). 




+250 KHz +265 KHz 



FIG 11-14 POWER SPECTRUM OF MIXER OPERATION 


04 
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The return power spectrum frequencies occur due to 
the forward viewing of the scatterometer This power is 
divided and inserted m the frequency spectrums about the 
two mixing frequencies. When the frequencies about 
250 KHz are subsequently sampled, this power division 
will cause an incorrect do calculation to occur. 

This method also eliminates utilization of those 
areas defined by the doppler shift of from 5-10 KHz. 

This occurs due to the doppler power impressed m fre- 
quency around 250 KHz and 265 KHz overlapping Back- 
scatter coefficient calculations will therefore be in- 
correct for areas defined by these data points. 

These effects nullify all advantages which might 
arise by utilizing this method of calibration. Thus an- 
other method of calibration needs to be implemented 

Instead of the present method of continuous cali- 
bration being used, calibration for each hardware com- 
ponent could be implemented This would complicate the 
system due to the necessity of adding extra hardware. 

Other methods of using' continuous calibration are 
available The 265 KHz signal could be fed directly in- 
to the circuit after 'the mixer This would eliminate the 
type of power split previously mentioned and allow for 
continuous calibration of all components the previous 
method excluding the mixer. 
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Probably the best method would be that of using an 

/ 

altered SSBM to achieve the type of operation shown m 
Fig. 11-15. Utilizing this method, the mixer output 
would be that of a spike at 265 KHz plus the doppler 
shift frequency around 250 KHz The power loss, ampli- 
fier flucuations, and system noise could then be ac- 
counted for and corrected The power split m the mixer 
would not occur 

Signal Re duction and Logic 

Once the signal has been received and the center 
frequency of the doppler frequency envelope has been re- 
duced, the signal is processed to obtain the backscatter 
coefficient (a 0 ) This processing or signal reduction 
section is shoira m Fig. 11-16. 

Signal reduction is required to implement the real 
time CRT display. This signal reduction provides the 
plus variance for each area viewed. 

The signal, consisting of the doppler information 
and the calibration frequency, is passed through three 
amplifiers which allow for amplification of 30- 70 db . 

The signal enters three circuits, one to allow for cali- 
bration, another to filter the desired doppler informa- 
tion, and the third to the flight recorder. 









265 KHz 


norizontally polarized signal reception circuit 



FIG. 11-16 SIGNAL REDUCTION CIRCUIT 
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The signal to the flight recorder enters a mixer 
which transmits only the lower side band thus bringing 
the center doppler frequency of an area to zero hertz. 
This is accomplished by utilizing a voltage controlled 
oscillator (VCO) to generate the frequency desired for 
the mixer’s local oscillation frequency 

The use of doppler frequency for area delineation 
necessitates a frequency shift for the electronic proces- 
sing A ground cell is defined by doppler frequency 
differences, for example 8064 Hz to 8300 Hz. This ground 
area would have a change m doppler frequency of Af^ = 

236 Hz (see Fig 11-17) Only the area information de- 
fined by the Af^ and the 3° beam width is desired. The 
VCO is biased to oscillate at the signal midpoint fre- 
quency of this area. 

f vco = 

= 2 50 KHz +-- ° t> g t83 ‘ l ° 

(2-i4) 

where f ( -f ^ = lines of constant doppler frequency 
; " shift which delineate a ground cell 

This frequency, f VCQ , for this area is fed into the 

mixer The mixer creates a signal of frequency spectrum 
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resembling the diagram m Fig 11-18 Part of the pre- 
vious frequency spectrum (i.e 250 KHz) is folded over 
onto itself causing the voltage amplitude m this region 
to add. 

The signal is amplified and passes through a vari- 
able low pass filter (VLPF) The amplifier is designed 
such that the information previously at 8182 Hz, now 
D-C, will be amplified, thus retaining information at 
this point This restriction can be eliminated if the 
information lost is quite small compared to the total 
information reflected from the area under consideration 
The D-C amplification is necessary due to the nature of 
the calibration techniques used. Another restriction is 
that the amplification must be almost constant over the 
frequency range from 0-1000 Hz for reasons to be shown 
later 

The frequency response of the VLPF is controlled 
from the synchronizer The frequency response varies to 
compensate for velocity and scan angle variations This 
VLPF variation is necessitated by the geometry of the 

CRT display If the frequency response of the VLPF was 

/ 

held constant, the ground area represented by the fil- 
tered frequency range Af^ would vary This occurs due 
to the nonuniformity of the lines of constant doppler 
frequency (see Fig. 11-19). 



Amp li tude 
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FIG. 11-17 

GROUND APPEARANCE OF LINES OF 
CONSTANT DOPPLER FREQUENCY 



FIG. 11-18 


MIXER SIGNAL OUTPUT 
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To keep the area constant for each cell at any an- 
gle, the frequency response (f ) of the filter (1 e. 
equal to -^-Af^ m Hz) must vary. To analyze only the fre- 
quency spectrum of the hypothetical area, the VLPF fre- 
quency response will be 

(y. - -J A-f^ =j(l3CHz) ~ I I? Hz (2-15) 

The response m only one -half of the doppler fre- 
quency representing the area since the frequency spectrum 
m f is the folded- over spectrum area from the mixer 

The information contained m f is equal to that Con- 
i' I 

tamed m Af^. The frequency response f^ is not likely 
to exceed 1000 Hz. Therefore, the amplification level 
of the amplifier m this circuit needs only to be flat 
m the frequency range 0-1000 Hz. 

The signal now enters a linear detector. This re- 
duces the signal to the voltage envelope as shown m Fig. 
11-20 

The reason for considering voltage here instead of 
power is explained through the operation of the calibra- 
tion system As mentioned previously, the signal from 
the three amplifier stages enters three circuits, one of 

i 

which is the calibration signal reduction section Since 
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flight direction 



FIG. 11-19 

SCATTEROMETER GROUND VIEWING 
FOR CONSTANT DOPPLER FREQUENCY 
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| received ~ 

transmitted 

f signal 
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l\ 

voltage 


Li 

L 

liUl/lA envelope 



FIG 11-20 


A LINEAR DETECTOR 
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the three circuits are m parallel, the voltage levels at 
all frequencies entering these circuits will he equiva- 
lent. Up to this point power conservation has been a 
lain goal. The voltage relationships now have to be con- 
sidered such that equivalent signals will pass through 
both the calibration signal and doppler signal reduction 
sections Loss of power, however, must be avoided. 

Power is now considered m one of its representative 
forms . 


P= VI = V-p 



(2-16) 


where P = power 

V = voltage 
P = resistance 
X = current 

If the impedance of the reduction circuits for the 
calibration signal and doppler signals are equivalent, 
a power ratio would be equivalent to 




2 


C2-17) 
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where d = implies doppler 

C = implies calibration 

This power ratio will occur later m the final op- 
eration of the calibration system. Therefore further op- 
erations will be concerned with voltage levels. 

The amplified return signal initially enters a mixer 
m the calibration signal reduction section. The cali- 
bration oscillation signal (265 KHz) is mixed with this 
signal. The mixer transmits only the lower side band 
The calibration signal m the return signal is therefore 
beat to zero hertz. The resultant spectrum is shown m 
Fig. 11-21 

The signal enters an amplifier equivalent to that 
m the doppler signal reduction circuit This amplifier 
must amplify D-C m order for the calibration signal to 
be subjected to the same treatment the doppler signal 
receives m its reduction circuit. 

The amplified signal passes through a variable low- 
pass filter comparable to that contained m the doppler 
signal reduction circuit The frequency response f of 
this VLPF is equal to the f m the other reduction cir- 
cuit. This causes the impedance values of both VLPF’s 
to be equivalent. 

The D-C calibration signal is now transmitted to a 
linear detector which also is comparable to the one con- 
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FIG 11-21 

RESULTANT SPECTRUM AFTER MIXER 


OPERATIONS IN CALIBRATION REDUCTION CIRCUIT 
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tamed xn the other reduction circuit. The return sig- 
nal is therefore subjected to equivalent operation, ex- 
cept for the beat frequencies m the mixers, m both re- 
duction circuits This means that the resistive effects 
of both circuits are equivalent and that the use of vol- 
tage vs. power handling is justified upon the division 
of the two resultant signals. 

The signals from the two reduction circuits enter 
an analog divider. Using voltage magnitudes, the signals 
will be passed through diodes whose output response will 
be the logarithmic value of the input. The two logarith- 
mic signals are then subtracted. This value is doubled 
to account for the voltage squaring m order to represent 
the power ratio The signal is then passed through an 
antilog device to obtain the final signal output 

There exists the possibility of using power magni- 
tudes for operations m the analog divider. Since the 
impedance values in both reduction circuits are equiva- 
lent, the signal power entering each of these circuits 
will be equal The power entering each of these two cir- 
cuits will be less than one-half of the previous total 
power due to power also entering the third parallel cir- 
cuit. After the signals are processed, the power con- 
tained m the calibration signal (a«P ) is an equal frac- 
tion of its -previous value as is the power contained m 
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the doppler frequency signal (a*P^) The analog divider 
must be designed such that the logarithmic values of the 
power magnitudes is obtained, subtracted, and the anti- 
logarithmic power magnitude acquired 


In a P - In a + In P 


therefore 

P 

In aPd “In Q = In Pd "In P c = (n— (2-i8) 

where Q = constant 

Whichever signal process is chosen, the calibration 
signal must be designed such that its logarithmic value 
represents signal variance m the system This requires 
that the signal magnitude of the reduced unaffected sig- 
nal will give a logarithmic value of zero This would 
indicated that the return signal was also unaffected by 
signal noise and losses If this final reduced calibra- 
tion signal is different from the specified value, its 
logarithmic value will not be zero Then, the division 

process will correct the doppler signal to its unaffected 

1 

value, assuming the doppler signal and calibration signal 


were equally affected by the system 
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After all gam variation has been removed by the 
analog divider, the signal enters the integrate -and-dump 
section. This section acts much like a charging capaci- 
tor in which many return signals are allowed to add. The 
signal will appear as shown m Fig. 11-22. When a cer- 
tain time limit is reached, the voltage value obtained is 
discharged and passed into the sample-and-hold section. 

The integration technique used here must actually 
find the average value of all the input signal voltages. 
The pulse rate frequency 45 KHz, or 225 KHz, determines 
the number of input signals. 

If the reception occurs for .2 seconds during the 
45 KHz PRF, nine thousand signal voltages will be placed 
into the integrator The voltage value sent to the sam- 
ple-and-hold section must be the average of all of these 
signal voltages. This is quite easily attained. An 
example of how this can occur is shown m Fig 11-23 
The signal voltages will charge C^. During the period 
of antenna lobe readjustment, no signal is coming into the 
integrator, thus, the voltage across capacitor will be 
accumulation of the 9000 signal voltages Then switch 
is thrown. The capacitor will discharge its voltage 
as 
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Vc,(t) = v mox e 


R.C 


where C = C- 


(2-19) 


The voltage will be allowed to discharge until 
exp(-^“) = g - QQQ then, the voltage across capacitor C 2 
will be one nine- thousandth of the maximum voltage and 
thus, the average voltage When this time occurs, 
switches and are thrown. Capacitor retains its'' 
average voltage value and capacitor C£ discharges rap- 
idly The process of capacitor recharging starts 
again. Prior to being thrown again $2 is closed, then 
opened, to discharge capacitor and thus prepare it for 
its next charge 

This is just one method of integrating- and-dumpmg 
plus sampling- and-holdmg Switches and S2 must be 

thrown m their proper sequence during the antenna lobe 
readjustment period This period is quite small and thus 
sets quite stringent specifications for this part of the 
system and its operations. Whatever method is used, it 
will follow these lines It must be noted that this is 
just one of the many possible methods of integrating and 
sampling 

The analog- to-digital converter has a switch m 
series with an extremely large resistance. The switch 
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1 2 3 4 Signal Reception 

FIG. 11-22 

VOLTAGE REPRESENTATION OF SIGNAL INTEGRATION 



FIG. 11-23 


EXAMPLE OF INTEGRATE AND DUMP CIRCUIT 
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will be thrown, thus sampling the voltage, without too 
much voltage dram. The signal is then m digital form 
for normalization purposes 

Using this average return power, the backs catter co- 
efficient can be determined by utilizing the radar equa- 
tion, [18] 


w s = 


1-1 2 if 


( 2 - 20 ) 


where V4 = averaged return power 

G 0 = antenna gam 

transmitted wavelength 
)= antenna pattern functions 
© - incidence viewing angle 
(j) = transverse viewing angle 
h = aircraft height 
T/ = aircraft velocity 
0^ = average backscatter coefficient 

This equation is rearranged to obtain an equation for the 

backscatter coefficient 


W s a^/r) 3 h 2 or 

p t G: r f,(e) ^(0,4*) 


( 2 - 21 ) 
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This operation occurs m the data normalizer section. A 
block diagram representing this section and its inputs 
is shown m Fig. 11-24. 

The received signal will appear as shown (Fig II- 
25) during the passage through the signal reduction sec- 
tion 

The final pieces of hardware m this section of the 
scatterometer are the synchronizer and its accompanying 
equipment This hardware constitutes the logic section 
and governs all scatterometer operations. A block dia- 
gram of this hardware appears as shown m Fig. 11-26. 

The synchronizer commands will determine the vari- 
ance of the determined Oo from its true value As shown 
previously, this variance is determined by the number of 
independent samples (N 1 ) This N.. is 

N s time to pass a resolution element of width Ap 
i time per independent sample 

where: time/s ample = - ^ r 

M Af d 

therefore 


^ Afd Ap 

c i/H V 


C2-22) 


where V = velocity of the aircraft 

The velocity and antenna patterns control the effective- 
ness of the system. Also, the increase of the reception 




signal 


FIG. 11-24 THE DATA NORMALIZER 








FIG. 11-25 SIGNAL VOLTAGE ENVELOPE SPECTRUM AF 



VARIOUS OPERATIONS 
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FIG. 11-26 SYNCHRONIZER INPUT AND OUTPUT 
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and integration times will increase N . This will also 
affect the effectiveness of the system. Thus, the design 
and operations of the synchronizer play a dominant role 
m the usefulness of the s catterometer 

Recorder 

The returning signal contains information from 
ground areas at all angles of incidence up to and m- 

i 

i 

eluding 50°. It is desirable that all of this informa- 
tion be retained. A tape recorder system shown m Fig 
11-27 is utilized to record on dual track tape informa- 
tion m analog form for subsequent processing One track 

1 

is open for the signal from the horizontal reception cir- 
cuit and the other for the vertical reception circuit. 

The signal is received from the reception portion 
and passed through fhe three amplifiers. A code oscilla- 
tion is injected into the third amplifier of the circuit 
that is receiving the cross polarization component The 
recorder is receiving either the return signal or the 
code oscillation. Thus, when the code signal is trans- 
mitted through the vertical reception circuit, the hori- 
zontally polarized return is being recorded When the 
code signal is transmitted through the horizontal recep- 
tion circuit, the vertically polarized return is being 
recorded. This is useful during post-flight evaluation 



horizontally polarized 
signal reception circuit 
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Circuit ’ r 


closed when 
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FIG. 11-27 RECORDER CIRCUIT 





62 


for piecing together the horizontal and vertical tape 
portions . 

The blanking circuits are used to eliminate passage 
of any signal during the period of antenna lobe adjust- 

i 

ments . This eliminates any spurious signal from the an- 
tenna, code oscillator, or system from being recorded as 
a signal during this period 

The line drivers are used for impedance matching 
purposes They are designed to operate over the fre- 
quency range of 150 KHz to 350 KHz with less than 1 % am- 
plitude and frequency distortion. 

This taped signal will be used for post-flight 
evaluation of ct 0 vs. 9 for those areas which this is 
possible This signal is m analog form on the dual tape 
track for purposes of distinguishing between horizontally 
and vertically polarized reception Post-flight evalua- 
tion will utilize the identical signal reduction and 
logic section of the scatterome ter The required adjust- 
ments will be made to allow for o 0 vs. 9 evaluation for 
9 other than 50°. The procedure for accomplishing this 
will follow from the information given about the signal 
reduction circuit. 
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CRT Disp lay 

After signal normalization occurs, a value of cro is 
obtained and transmitted to the cathode-ray- tube display 
circuit (Fig 11-28) 

The display circuit receives the signal as a vol- 
tage representation of Oo The information from all the 
areas on a single sweep appear on the display as the top 
line after having been released by the scan storage and 
passed through the digital- to-analog converter. The ana- 
log signal is necessary for continual CRT sweep This 
information was also transmitted back into the gate which 
will then release information from the newest sweep and 
subsequently this returned information. The analog in- 
formation is used to intensity modulate the signal of the 
CRT display As the information is transmitted to the 
gate and released prior to newer information, the dis- 
play shows the a 0 information declining, hence the fall- 
ing raster technique. This will appear as shown m 
Fig. 11-29. 

Since a great deal of effort is expended m evalu- 
ating Co as correctly as possible, the screen brightness 
should adequately represent the signal voltage value. A 
curve of input voltage vs. output brightness m lumens 
will appear as shown m Fig. 11-30 
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FALLING RASTER TECHNIQUE 
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FIG. 11-30 

BRIGHTNESS VS VOLTAGE 
CURVE FOR CRT DISPLAY 
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This curve m non- linear. Therefore a gamma- correc- 
tion factor is required in the display circuitry to ad- 
just the input voltage. This causes the curve to m ef- 
fect become linear, thus adequately representing o 0 
through brightness [19]. 
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CHAPTER III 

MEAN DETERMINANCE OF A BACKSCATTER 
SIGNAL USING A HALF-WAVE LINEAR DETECTOR 

Signal backscatter is subject to fading. The signal 
reflected from the ground and received at a point above 

the ground will have amplitude variations commensurate 

/ 

with slight variations of the receiving point This oc- 
curs because the signal received is the vector sum of all 
the individual returning signals existing at that point 
Thus, at a different point, the vectors will be different 
and will produce a different sum. Therefore, for an air- 
craft flying over the reflecting ground, the return sig- 
nal will be subject to fading. 

The s catterome ter described m Chapter II integrates 
the return signal, from a ground cell having a doppler 
frequency difference of 280 Hz viewed at an angle of 50°, 
for periods of 0 2 seconds The scatterometer will inte- 
grate and average numerous transmitted, and subsequently 
reflected, pulses. If the pulses reflected from a spe- 
cific ground area were averaged over an infinite length 

\ 

of time, a true mean or average signal value could be 
found. Since infinite time averaging is impossible, it 
becomes necessary to find how the finite time average 
approximates the true mean signal value reflected from 
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the ground This involves determining the variance or 
decibel range about the true mean signal value within 
which the finite time average is located. These values 
will be found for the operations of the s catterome ter 
described m Chapter II. Also, possible improvements for 
the s catterome ter averaging are given. 

Fading 

If the transmitted signal wavelength is large m 
comparison to the ground distance required for appre- 
ciable slope change, the ground is considered "smooth". 
The signal reflection from the smooth surface will fol- 
low Snell’s Law of reflection. Snell’s Law depicted in 
Fig. III-l, states that a ray incident upon a surface 
makes an angle with the surface normal which is identi- 
cal with the angle the reflected ray makes with the sur- 
face normal. This law is the basis of all geometrical 
optics observations It applies only for smooth sur- 
faces. The determination of a surface being smooth de- 
pends upon the individuals definition of large and appre- 

» 

ciable m reference respectively to wavelength and slope 
change. 

Similarly, if the signal wavelength is small m com- 
parison to ground variations, the ground is considered 
"rough” Statistical distribution methods are used to 




angle of incidence , 9j, equals angle of 
reflection, 0 2 , for a smooth surface 


FIG. III-l SNELL’S LAW OF REFLECTION 
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determine the signal reflection or scattering from a 
rough surface. These methods will define the surface 
roughness and will determine the distribution of the am- 
plitude and phase of the scattered field. Nearly all 
natural terrain surfaces can be assumed "rough”. 

Statistical methods, which define a rough surface, 
describe the distribution of the surface height devia- 
tions from the mean or average level with a density curve. 
One type of density curve, the Gaussian, often arises as 
the solution to scattering from rough surfaces The 
Gaussian density function is described for a one dimen- 
sional surface as 





(x-m 2 

7F 


(3-1) 


where V)** the mean or average height 
f(X]= the density function 
XT 2 = the variance about the mean 

As the variance, o 2 , approaches zero, the Gaussian 

distribution will bunch about the mean value until, when 


.2 _ 


= 0, the density function will be a unit impulse lo- 


cated at the mean value. 

Another density function which is used to describe 
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surface and signal amplitude distributions is the Ray- 
leigh density function, given by [20] 



C3-2) 


Some properties of these curves are that they are 
non-negative and that the area under the curve is always 
equal to one. 


oo 


fcxjolx - I 


-CO 


The mean value ri can be found from the density curves 
using the expected value technique, where E{ } is the 
expected value 


<*o 

E = I x 

£ £ X n j = j X^C*) dx 



The variance, or amount of deviation around the mean 
value, is used to define the range of values that the 
distribution assumes. It is defined as a 2 . 
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& E [tf- = E{tf - E z {x] 


Signal fading necessitates time averaging the sig- 
nal m order to evaluate the mean value of power re- 
flected from a surface type. This mean value is uti- 
lized to provide an estimation of the overall roughness 
of the ground area viewed. Knowing the variance that 
the averaged power value has from the mean value, the 
approximate backs cat ter coefficient can be calculated 
[ 21 ]. 



f, (e) f z (e,<t>) P t l 3 


(3-3) 


where 

( X = backscatter coefficient 
averaged return power 
h - aircraft height 
If- relative velocity 

Since all other values can bi 

length of time, the variance 


= transmitted power 
G 0 = antenna gam 
~X = signal wavelength 
■f C ) = antenna paterns 

held constant for a given 

of the calculated signal 
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mean from the true mean will be equivalent to that of the 
averaged signal value plus noise, thus, 

£r 0 = + n 


where = variance of the obtained backs catter coef- 

ficient from its true value 

£ cJj = variance of the averaged b acks cat te red 
J power from its true value 

n = averaged noise 

The averaged noise will be assumed to be small com- 
pared with the signal, therefore, 

~ cT 0J S 

The return signal will be the result of transmitted 
pulses thus giving the effect of the return signal being 
sampled at the transmitted pulse rate. The Central Limit 
Theorem implies that the amplitude distribution of a sam- 
pled signal being averaged can be considered Gaussian. 
Therefore, the variance of the averaged power from its 
true mean will be determined using statistics involving 
Gaussian distributions 
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Detection 

Consider the case where the return signal is mixed 
with another frequency, passed through a low pass filter, 
and detected with a linear detector, as shown m Fig. 

Ill- 2. The spectrum of the signal will change as shown 
m Fig. III-3. 

The spectral density of the output of the linear 
detector, S^ff), is evaluated from the input spectral 
density S (f) , where 



w he n | f | 5 A 
O'f'Inc.f'Ml 5 e. 


(3-4) 


where 


A - 
Afd “ 

m- 


amplitude 

doppler frequency used to define ground 
area 

impulse function 


and [22] 


s,m= 


s. 

Znr 


fff) 


— oo 


(3-S) 
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FIG. Ill- 2 RETURN SIGNAL CIRCUIT 



FIG. Ill- 3 


SPECTRAL CHANGE IN RETURN SIGNAL CIRCUIT 
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The variance, cr^ 2 , of S^ff) is the area under its spec- 
tral density curve. 



Substituting this value for a 2 , 

X 


S/f) 


5y(f)- 


27 r 


ttf) + 


i 


becomes 

jVi S/f-f'Xf 

-0O 


+ 


T 3, ff) 


(3-7) 


The last term of the output spectral density equa- 
tion is the convolution of the input spectral density. 
Thus, 

OO 

J5 x (fjs x (f-r;df' 

— CO 



IF 1— ■A'Q 


Substituting this into eq. (3-7), S (f) is finally ar- 
rived at m terms of frequency, f, and amplitude, A. 
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Slf) = 


~h 


ushe-n lfl~'/z A ^ 

O o^io einjise. 


+ 


O 


!j£L 

Fet J wlaet' 1 ifl -^^4 

o -^her-Wise 


(fl &fdO &£,) wkenlfl-A^ 


(3-8) 


The spectral density after linear detection, S^(f), 
will appear as shown m Fig. III-4. 

The autocorrelation function, R (T) , which is used 
to determine the amount of independence between two parts 
of the detected signal, is obtained by taking the Fourier 
cosine transform of the spectral density S (f) [23] . 

oo 

Ry(T)= ^ cos(ZftfT) d f (3-9) 

_ oO 


where 'Y = a time difference variable 

By substituting eq. (3-8) into eq. (3-9) the auto- 
correlation function becomes 
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/ V 7> "^r + 2 j ^/Co 3 (2^fT)df 

O 

+ 2 5^('-^) co5( ^ fT)df 

(3-10) 

Integrating eq (3-10), using integration by parts, 


R (T) is derived as 

y 


R,m- 


A Afd 

27T 


-f- 


5/^1 T 

2 2/rT 


_£} l!j2 T 

ZK ZKT 


a£| 

ZXA-fj c^C2/rf Odf 


_ $ a /3 s j v ft a fd T 

2K 2 S/tT 

-f A - A cos 2 tt a-Q 7" 

(2?r) 3 T* 


(3-11) 


The trigometric identity for sin 2 9 is used to sim- 
plify R y (T) 

2 sm z 0=1- cos 2 0 


(3-12) 
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The final result for the autocorrelation function 
of the signal after linear detection is therefore 


Q (T)- — /-] Slfl T R SlfflXkfj T 

(3-13) 


The first term of Ry(T) is the square of the mean 
value of y(t) (i.e. riy 2 ) The rest of Ry(T) evaluated 
at T = 0 gives the variance, Oy 2 . 




R 


8 ± £a . JH} | + _l ) 

H n ^ ' K) 


(3-14) 


This follows from the fact that Ry(0) is the mean square 
value of the signal, y(t), transmitted from the linear 
detector. This value of R (T) is obtained by assuming 
that the mean value of the input signal X(t) is zero, and 
that X(t) is a sample function of a real Gaussian random 
process . 

In light of the criteria concerning the return sig- 
nal, this assumption may at first seem erroneous There- 
fore, an insight into the statistics of the return signal 
will be useful The amplitude of the return signal is 
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Rayleigh distributed. Let the signal be interpreted as 
the sampled return of a continuous wave CCW) system in- 
stead of the actual return of the transmitted interrupted 
continuous wave (ICW) signal. Then, the amplitude and 
samples of the signal would appear as shown m Fig. III- 
5. 

The amplitude distribution of the sampled signal is 
Rayleigh distributed. However, the distribution of the 
mean value obtained by averaging M independent samples 
will tend to become Gaussian as the number of independent 
samples N increases. 

d, + QL^ +■ • * * + O-M 

N 

where the sampled mean 

Ct<_ = the (.th voltage sample 
N ~ the total number of samples 

As N increases, the sampled mean will approach the true 

mean, n, of the signal, and the probability distribution 

of the sampled mean, f. (x) , tends to be Gaussian re- 

a N 

gardless of the probability distribution of each inde- 
pendent sample This is a statement of the Central Limit 
Theorem which is valid for independent samples having a 
finite mean, n , and a finite variance, o 2 [24] 




time 


signal before 
linear detection 


FIG III-5 


time 


signal after 
linear detection 


AMPLITUDES 
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Gaussian techniques are valid for determining the 
autocorrelation R^(T) of the return signal pulses. It 
must be noted though, that the probability density of 
the sampled mean is not truly Gaussian unless N ap- 
proaches infinity. For all practical purposes, consid- 
ering a Rayleigh distributed input amplitude, the sam- 
pled mean can be considered Gaussian after about ten in- 
dependent samples [25] 

In using Gaussian techniques, trouble arises due to 
the mean value of the incoming signal, n , not being 
identically zero, an assumption necessary m determining 
the equation for S (f) as a function of S x (f). This can 
be corrected for by placing a stringent specification on 
the calculated R (T) as to what constitutes an indepen- 
dent sample The qualifications for sample independence 
are arbitrary but they must be strong enough to account 
for any error which will arise due to the inherent in- 
correct assumption made m order to calculate R^CT). 

Before these qualifications are stated, the auto- 
correlation function will be converted to a more con- 
venient form A new function, K (T) , is formed by sub- 
tracting the constant term from R^CT). 


K 5 (T)= R^CT) - 


ZK 


(3-15) 
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Substituting for R (T) , K (T) becomes 


K 3 (r)= 


f) sin T 




/? Sf^ 2, K A-Pd T 
{2K) Z K Afol T 2, 


C3-16) 


To rid K (T) of the amplitude value A, the function 
K (T) is divided by its initial value, K^.(T = 0) , thus 
obtaining a normalized value p(T) 


f(T) = 


K s tO) 


K„ ( T) 

V 0) 

Pt_ A -Pol ( I + 7f ) 


The equation p(T) finally becomes 


C 3 - 1 7) 


sin k A-& T sm z K Afj 

KAfjTU*#) (znr)* T 1 (ff+i) (3-i8) 

In some notation, fy (T) is called the autocovari- 
ance function, and p(T) is called the normalized auto- 
covariance function, which has also been called the auto- 




correlation function [26] 

The autocorrelation p(T) goes to zero at 
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A-f<( T= I, Z, 3, ,n, Therefore X n - 

A good criteria for independence is for the auto- 
correlation function to be less than one-tenth after a 
certain time T. For convenience, the criteria will be 
set as 


£(T) - ft 2 for. all T > T, 


and 


p(T 2 )=0 wkere % > T t 


Using these conditions, 



sin(nff) sm z (nn) 
n(fr-i-i) + (2^) a n s (ft+l) 


C3-19) 


p(T) /^■a when ^ T n 


and 
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The criteria n = 3 will satisfy these two conditions. 


r- 3% - 1. 


where tZ/ 5 = time between independent samples 

Under this criteria, sampling at every seconds 

A *4 

will give a set of independent samples. If the sampling 

time is slightly greater, two adjacent samples will have 

some correlation which is less than — , which also satis- 

K 1 - 

fies the independence criteria If the sampling time 

3 

is less than — r seconds, only the string of samples be- 

^ 3 

mg equal to or greater than seconds apart can be 
considered independent 

The first sample is the first independent sample. 

3 

Therefore, at seconds later, the second independent 
At/ 

sample can exist, for (L r 2^ the third, and so on 


Sampled Variance of the Mean 


The amplitude distribution of each pulse received 
is random m a Rayleigh manner. The Rayleigh distribu- 
tion has the mean and variance shown. 


fcx)=-~ e ^ Ucx) 


( 3 - 21 ) 
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OCJ p 

Yj - X fcx) d x = \ 


* i 

iL p 2 

o<2- ^ 


(3-22) 


where c>c = value where f(x) is maximum. 
U(X;= unit step function 


variance 


QO 

<o~ 2 = 5 x£ - o : 


= ( 2- t)°< 


(3-23) 


For independent samples, a , the sampled mean , will be 


An “ 


Ct, 4~ Ci* ■b • • • +■ CLf 


For independent samples the variance, cr 2 ^, is 

GT* = C^ z 4- cr 2 + ... + rr i 

0..+ Cl,+ •• a i a, Yet.. 
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For any signal from the same patch of ground, every 
return pulse will have the same variance with identical 
Rayleigh distributions. Thus, 





L 


and 

CC+ fe** ■•■tC _ A/0~ z = l T z 

J\J Z n l n 

(3-24) 

With an increasing number of independent samples, the 
variance of the sampled mean goes to zero, i.e. the sam- 
pled mean approaches the true or Rayleigh distribution 
mean 

From Tchebycheff ' s Inequality, the probability P 
that is within a certain range is given by [27] 



< f) N - - 



2 . 


N c 2 - 


(3-25) 
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Since the incoming signal, x(t) is Rayleigh dis- 
tributed, the ratio of the mean squared, r\ 2 , to vari- 

Jv 

ance, cr 2 , is 

A 

t £ 


.272 

IT 



Thus , 


.272 

N c z 


C3-27) 


where 


2 / - K 



.£72 

yc 2 


Consider finding the probability that the value 
is within ± 3db of ri x where A^ is the averaged value of 
the return envelope voltage V . 


3 c lb = 20 locj 


A /). 





a 4= 7 X C'«) 


(3-2 8) 
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This requires that 


A ft N 

% 


0*0 >?, 

% 


- ( i fc) -7.^/3 


C3-29) 


Therefore, solving for c, 

C =.03 

// _ jUA - h£l 

N (.H 1 3f N (3-30) 

To be greater than 99% certain that is within 3db of 

1 59 

h x , i.e. 1-K = 99, K = .01 = — , therefore N = 159 

Thus, 159 independent samples are needed Consider, the 
probability that A N is within ±ldb of q . For this, 


Idh = 20 log 

= 2 0 log ( I +■ c) 


Thus , 


i*c = 1.121 

C = .! 2 I 


(3-31) 
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Therefore 

.272 . ! 2 .b 

K " N(-I 2 I ) 2 " N 


To be greater than 90% certain that is within ±ldb of 
n , using eq. (3-27) and eq (3-31) 

A 

I - K = - c ?0 => K = -10 

K-.IO-- 2 — =^>/?& = N 

To be greater than 99% certain that is within ±ldb of 
”x 

I- K = .1? =5> K = .Ol 

t 

K=.0/ = -^p- => N= ISloO 

Consider scanning a cell area having a change m 
doppler frequency of at least 280 Hz. An independent 
sample will occur m the time specified by eq. (3-20) 

3 3 

t, s = ^ = geo r ' °107I seconds 
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In viewing the same cell area for 0.2 seconds, which 
is the cell viewing time for the scatterome ter system 
described m Chapter II, the number of independent sam- 
ples will be as shown 



T_ 

t ls 



■ ,2 

.01071 



1 1 


(3-32) 


Thus, 19 independent samples are received and the prob- 
ability shown m Table III-l is derived with N - 19 from 
the probability equation, eq (3-27) 




■ 272 1 oH3 

C l N C* 


Table III-l shows that for 1000 values of V s ob- 
tained from viewing for 0.2 seconds, cells having a 
change in doppler frequency of 280 Hz, at least 995 of 
the values are correct to within ±8 6 db of the actual or 
true value, rj , at least 990 of them are correct within 
±6.82 db, at least 980 are correct within ±5 34 db , etc 



The Tchebycheff Inequality- 
Using 19 Independent Samples 
and Rayleigh Distributed Signal 


Percent of 
Certainty 

A^ Within 
±K db of ri 

Value of C 

0 0% 

0.98 

. 1195 

10% 

1.04 

.126 

90% 

2.78 

. 378 

95% 

3. 72 

.535 

96% 

4.04 

.59 8 

9 7% 

4 56 

.690 

9 8% 

5.34 

. 845 

99% 

6 82 

1 195 

99.5% 

8.60 

1.691 


TABLE III-l PROBABILITY LIST 
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Since, for the conditions stated, no averaged value can 
be said to be correct to within less than ±0.98 db of the 
true value, this constitutes the lower limit for accuracy 
The number of independent samples is given by 


N 


T 
t. . 


+ 



(3-33) 


where ~f~ = viewing time 

doppler frequency change across a cell 
"t Li ~ one independent sample time period 

Therefore, either an increase m T or an increase m Af^ 

would improve the mean determmance of V s for this system 

The above equation for N holds assuming the criteria that 

sampling occurs once every t seconds. For sampling 

faster than this, equation (3-33) reduces to 



X 

t,s 



(3-34) 


Due to this faster sampling, one independent sample will 
exist for every period duration of t g seconds. Thus, 
m averaging, numerous pulses will represent one inde- 
pendent sample 

There exists another method for finding the prob- 
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ability, P, that a value z lies within a certain region 
(a,b) on the density curve f (x) , is given'by [28] 

b 

Pjai2it,]= 5 • 

CL 

Using the Central Limit Theorem, the distribution 
of the sampled mean, A N , about the mean, n x , of the m~ 

I 

coming signal, x(t), will be Gaussian Thus, using the 
definition of probability, the probability that A^ will 
be within a certain db range of the mean will be 



The db range is given by 


dfc>= -ZO log ±20 loj(n-c) 


(3-36) 
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The probability integral is now reduced to 





where 


(3-37) 







W 


This is further reduced using the quantity called 
the error function, erf [29]. 


erf x 


/ 

VsF 1 




e 2 dz 


(3-38) 


The probability of being within a certain range 
is now given by 


P = 2 erf Jb = i? erf 


(3-39) 
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The db range for this probability can readily be 

evaluated from the above equation The variance, a 2 . , 

A N 

and mean, ri , values are the same as given before 

X 



± = fC 

2 . 



N 



oC 




(3-40] 


Therefore the probability will be 


P - 2 erf cVn 1-313 


( 3-413 


For to be within ±2db of rj^., the requirement 
is that C - . 260 


P= 2 erf (,240)(l. ^ I3)VF 


( 3-423 


when N = 4, P = 68% 

N = 9, P = 86% 

N = 18, P = 96% 

N = 25, P = 98% 

N = 36, P = 99.3% 

For A^ to be within ±ldb of n x > the requirement is 
that C = 121 
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P- 2 . Crf (- / £ ) X I* 13 ) VTT' C3-43) 

when N = 4, P = 35% 

N = 9, P * 50% 

N = 18, P = 67% 

N = 25 , P = 75% 

N = 36 , P = 831 

Determi nance of the Mean Power Values 

The previous calculations have determined the vari- 
ance of the sampled and time averaged voltage mean value 
from the true envelope voltage mean value, n . In de- 
termimng a value such as the backscatter coefficient, a 
value for the mean power is required Power is a func- 
tion of the envelope voltage squared. 

V 2 


where V ~ voltage 
P = power 

P = resistance - this value is made to equal one 
A mean power value W s can be derived by squaring the 




voltage mean value A^. 
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V\/ 5 =( fl N f 2L -ff ( 3 - 44 ) 

c= I 

where Cl t = t ^ ie tth independent voltage sample 
Another way to obtain a mean power value is by averaging 
the power values. 



2 


I 



( 3 - 45 ) 


where it can be shown that 

b n ^ K ) 2 

As N increases, confidence B N -> confidence (A^) 2 , 
therefore, the statistics of will be used to evaluate 
the statistics of (A^) 2 . This will give approximate 
values only. The statistics of (A ^) 2 are given after the 
statistics of B^ 

R = X 2 


Let 
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where X = Q-l = independent voltage sample 

The density of X, f (x) , is Rayleigh and given by 


-f„(x) = 


y _ ^ - 

_ ^ 2°< s 


o<‘ 


U(*) 


vVF * <r*‘=(2-?w 


(3-46") 


The density of R, f R (y) , is found as shown [30]. 


f ^-77n W) + 


(3-47) 


Substituting eq. (3-46), f R becomes 

i rvy - -i. 

it 


' 


+ 


-Vm 7 


2 e ** Ucvt) 




£ U("^T) U(y) (3-48) 


but 


U(-VTp --0 


since the unit step U(c) is defined as zero for negative 
c. Therefore, 
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(3-49) 


The mean and variance of R can now be found 


mean 


o o 


V E hl \ 5 'i ^ 


z <=< 


' jye" to “d 3 





E [f] = d a 

co y 



- ? °< H (3-50) 

variance 

<r/- E[f]~E'h} 

- 8«" - (£«*/ 

- ^ C<^ 


(3-51) 
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where = mean value of R 

(j"^= variance of R 

Using the Central Limit Theorem, the density func- 
tion, f R (z), of the average B w becomes Gaussian as N m- 
creases 

e 

(3-52) 

t 

where 



Thus, the level of confidence or percent of prob- 
ability that B n is within a certain db range of the true 
mean is, m the manner it was for in eq. (3-39), given 
by 

2erf 



(3-53) 
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Using the statistics derived for B^, it follows that 




2j*! ~ 

VaT 



and 


d b - -10 locj U + c) 


(3-54) 


For 
that C = 


B^ to be within ±ldb of q , the requirement is 
.260, giving the probability shown 


p = 2 erf (. 2 (oO)^ffT 


when N = 2, P = 40% 

N = 9, P = 56% 

N = 18, P = 73% 

N = 25, P = 80% 

N = 36, P = 88% 

N = 49, P = 93% 

The statistics of (A^) 2 are summarized below m the 
same manner that those of B^ were presented Now, from 


eq. (3-44), 
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N 





N 


By the Central Limit Theorem, the density A^, 
Gaussian 


£ a Cx), 

a n 


is 



l 

-V2?r ‘ 



c 



(3-56) 


Let 

5 = C Qu? 

VI 

, c r 

now, the density, f (y) , o£ S is, as m eq. (3-47), 
given by 





(3-57) 
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Substituting eq (3-56) into this equation, f g (y) 


becomes 


f s (j> = 




-UT-nf 

e 

+ e a<ri 


Ucj) 


(3-58) 


Squaring the exponent and separating parts, f g (y) becomes 




-tLiOil 

ITT 2 ' 


M -K 1 
0^+6 


-4 


- e 


zcr 


/5Fi r 2 - 1 LYT 


1 


COS] 


U^)j 


(3-59) 


(3-60) 


The mean value, p , of f (y) is given by 




% = J 3 ^ )oK f 


- oo 


(3-61) 


Substituting m £ g (y) , n s now becomes 
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It izw* 


CO 



O 



z$(sr*+fr)e 



= ■y* + -f 

= ^ (“Tp + l) (3-62) 


The mean value of (A ^) 2 , r) s , finally becomes 



(3-63) 


Finding the percent of certainty or level of con- 
fidence that (A ^) 2 is within +K db of n s is accomplished 
by performing the following calculations The db range 
m terms of c is given by 

K db = ± 10 log (1+-C) 


(3-64) 
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The probability P is, as shown m eq. (3-35) 


P{ %U+C)j 

- { (3-65) 


Substituting f (y) m eq. (3-58), P becomes 


P-- 




u, + o 

^ J 


/wf j vy 

%0-c) 

+ e 


r - tv ?-'?). 1 

e 


ZT' 


?sr L 


d 3 


(3-66) 


x 2 and dy = 2xdx. 


__ tx+nf 

Q d?C (3-67) 


now, let 



for integration purposes, let y 
There fore , 


P- 






e t 


4% 


l-C) 


dx = di 
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and 


VI = dx = 0" dw 


The level of confidence now becomes 

,z 


? = 


I 


o 1 +/ i 


-fzn 


e 


jr 

z 




d^ 


+ 


1 


V _ Ja7 

e £ diA; 


~ (3-68) 

Using the error function definition of eq. (3-38), 

P becomes 


p= erp^f^-) - erf 


+ er 




V'fcOd -■? 


V 


+ er 


-JWV-i 


tr 


(3-69) 


Substituting the value for ri m eq (3-62), P 


becomes 
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For (A^) 2 to be within ±ldb of n s > C = 
found from eq. (3-64). Let N = 36 



260 as 


JL = (3-71) 

r 


The probability that (A N ) 2 is within ±ldb of p s for 
N = 36 is 


P= erf(z 2 .l)-er-f( 2 n) + erf(i.3?) 

- <8 rf (" I- (o I ) (3-72) 

For large 6, i.e. 8 > 5, 

erf (3 erf (-p)- -erf ^ 
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therefore erf (22.4) = erf (21.4) and the probability- 
equation reduces to 


P - erf ( i. 3q) + erf(l.tl) 


= 21 , 7 . 


(3-73) 


which is approximately the probability for when 
C = .260 and N = 36. 

Now consider the probability that (A^) 2 is within 

±ldb of n . For N = 49 and C = .260, the ratio 5- becomes 
s o 

71 - j h H z ' = I 3 q 

q- y 4- ft- IJ- 1 

— = 5\ (via X I O’ 3 (3-74) 

T 

Therefore the probability becomes, with erf (a > 5) = 
erf (b > 5) 

P = erf ( I. LZ) + erf ( l.?8) 


II 7 


(3-75) 
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Summary 

The results obtained apply to the scanning scat- 
terometer described in Chapter II. The s catterome ter re- 
ception of the Rayleigh distributed voltage return signal 
is shown to contain 18 independent samples during the 
0 2 second integration time for a doppler bandwidth of 
280 Hz The variance of the backs catter coefficient from 
the true value was determined m terms of db variance and 
level of confidence, and their dependence on integration 
time was calculated A comparison of the average vol- 
tage squared, (A^) 2 , and the average power level, B^, 
shows that these two quantities do not tend towards the 
same values even though their confidence levels are 
nearly identical. 

The results obtained by the s catterome ter , aver- 
aging only 18 independent samples, i e. 0 2 seconds 
and 280 Hz, are shown to be inadequate due to the low 
level of confidence. Adequate results can be obtained 
by increasing the number of independent samples to 49 
This is accomplished by either increasing Afp the dop- 
pler bandwidth, or the viewing time for a cell. Doing 
such would cause improvements m the s catterome ter ac- 
curacy for determining the scattering coefficient, but 
the results of doing so would decrease the resolution. 
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CHAPTER IV 
FREQUENCY AVERAGING 


Frequency Agility 

Most radar s catterometers utilize monochromatic 
transmission. There exists the possibility that fre- 
quency agility or panchromaticity may mhance a systems 
ability to determine the scattering coefficient 

A frequency agile system works on the premise that 
there exists a frequency change, Af , which, if imple- 
mented between pulses, would constitute independence be- 
tween the two pulses. Thus, the transmitting frequency 
would be increased Af c between each pulse until the de- 
sired number of independent samples is reached However, 
there exists difficulties m using frequency agility 
The main problem can be seen by viewing the standard 
radar equation used m the scatterometer data normalize! 
to calculate the backscatter coefficient 


Ws ziVrf h 2 ^ 

P t G :f(e)f 4 (e» V 


( 4 - 1 ) 
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The variation in cro with respect to wavelength, X, 
can vary from X +2 to X” 6 dependent upon surface condi- 
tions [31]. For purposes of calculations, a o 0 depen- 
dency of X” 3 will be assumed If m using frequency 
agility, the wavelength range difference would be six 
percent of the center frequency wavelength (i.e. AX = 

. 06X), the variation m o 0 would be 


i 6<r a = i 


t 


( I + 


a 

K 


) 


= .0 8 S' 


(4-2) 


or 


£^7 ' - . 36 olb (4-3) 

where & A = variance m center wavelength = ±0.03 Xo 
= center wavelength 

SO^ = variance m backscatter coefficient 
Thus, a scatterometer having a transmitting center 
frequency of 13 3 GHz would have an inherent error of 
±0.36 db with a frequency change of ±400 MHz. Whether 
or not such a large frequency variance is necessary will 
be shown later. 
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The Autocorrelation Function Vs . Frequency Agility 

The autocorrelation for a low-pass spectral input to 
a linear detector has previously been calculated. This 
autocorrelation function holds for correlation between 
pulses having identical transmission frequencies. Pulses 
received having had different transmission frequencies 
will have a different autocorrelation function to deter- 
mine independence between pulses of different frequencies. 

Suppose that there exists a M pulse repetitive se- 
quence m a frequency agile scatterometer , and that each 
succeeding pulse frequency is increased by Af c , such that 
each returning pulse constitutes an independent sample. 
Then, M independent samples would exist every t sec- 

1 • S • 

onds . The interval t is dependent upon the mono- 
chromatic autocorrelation function, and it designates 
time difference between pulses of identical frequencies 
necessary for those pulses to be considered independent. 
When the time necessary for the M pulse sequence is less 
than t seconds, only one pulse at each specific fre- 
quency within each time period of seconds will 

constitute an independent sample (see Fig. IV-1). Thus, 
when three pulses exist m the M pulse sequency, l e. 

M = 3, three independent samples exist every t g 
onds 


sec- 
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The frequency autocorrelation function of a back- 
scattered signal has been shown to be [32] 

Z7C L sm Q l 

Sm\ c 

L sm 9 a£ 

where C - speed of light 

frequency autocorrelation function 

L = length of object being viewed 

© - angle of signal incident upon the object 

Af - frequency difference from one pulse to 
another 

Independence between signals of different frequen- 
cies will be assumed to occur when the frequency auto- 
correlation function is less than one- tenth. For con- 
vience of calculations, independence will be defined as 
occurring when p(Af) £ Therefore, letting the de- 

nominator equal 7T 2 , the calculation of Af c yields 

ZK L 5inO A~£ 
c 

, r - c _ 150 MHz 

c 11 1 s\n © l sino 





(4-5) 
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The ground length viewed, L, can be estimated using the 
diagram of Fig IV- 2. 

The doppler frequency value along the ground m the 
flight direction is given by 


I -_i^r 5m0 - sL y £ 5 1 n 0 

A C 


( 4 - 6 ) 


where +4 = doppler frequency 
if = aircraft velocity 
•p = transmitted frequency 
C = speed of light 

Q = angle of signal incident upon L 
Holding velocity and frequency constant, the deriva- 
tive of f^ will be 


df d = ^£cos 8 ole 


( 4 - 7 ) 


For small d6, this equation can be approximated as 

Af d = - cosQ aQ ( 4 - 8 ) 


where* Afd = doppler frequency shift contained along 
the distance L 

A0 = differential angle through which L is 
viewed 
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N = number of independent samples 


FIG IV- 1 

EFFECTS OF FREQUENCY AGILITY UPON SAMPLE INDEPENDENCE 



FIG. IV- 2 SIGNAL TRANSMISSION PATH DIAGRAM 
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The radial distance, r, from 0 to A (see Fig. IV- 2) 
is found by the relationship of a right triangle. 

r s _b. 

1 CO 5 0 (4-9) 

where h = height above ground for aircraft 

'f = radial distance from aircraft to ground 

For a small A6 , the arc AB can be considered a 

straight line. The arc has a length of AB = rA0. When 

the arc is considered as a straight line m the right 

triangle ABC, a relationship exists between L and A0. 


A B = r A© = L cos© 


(4-10) 


Substituting the relationship derived for r, the 
equation for A0 becomes 


A 0 = 


,L cos © 

r 


L cos z Q 

h 


(4-11) 


The equation derived for the doppler frequency shift 
contained along the distance L was 



iAT-f cos & A© 

c 


(4-12) 
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Substituting A0, the doppler shift equation becomes 

Afd tos 3 e (4-i3 


It is necessary to have an equation of L as a func- 
tion of other parameters This is obtained by rearrang- 
ing the doppler frequency shift equation to obtain the 
following 


h c 

2 v -f cos 3 0 


(4-14) 


where all parameters are as previously defined 

A value for L is necessary to solve for the fre- 
quency change, Af c> necessary to afford signal indepen- 
dence between signals of different frequencies. Sub- 
stituting L into this equation will constitute the final 
step m the derivation of the Af c equation Thus, 



c 

2 L sin Q 


v cos 3 Q 
h Afd Sin G 


(4-15) 


This is only an approximation which tends to hold 


as A0 becomes small. 
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Consider the case where extreme conditions for a 

scatterometer viewing from an aircraft are 

maximum aircraft velocity, v = 150 m/sec 
minimum aircraft height, h = 600 m 
minimum Af^ = 280 Hz 
0 = 50° 

frequency, f = 13 3 GHz 

The necessary frequency change, Af c , needed to af- 
ford signal independence becomes, with substitution of 
the given parameter values 

r (I5-0)UH25) 

Afc= (465)(2io)(:7uJ ,i - 3 10 

= q.n * /o 4 

= q.i i MHi 


Thus, for the extreme conditions stated, a frequency 
change of 5 MHz implemented between pulse to pulse trans- 
mission will assure that the reflected pulses will be 
independent . 

Previous work has centered on applying all tech- 
niques to a scatterometer having a transmission frequency 
of 13.3 GHz For this scatterometer, transmission for 
0 2 seconds and reception from a ground area delineated 
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by a doppler frequency change of 2 80 Hz m the flight 
direction caused 18 independent samples to be received 
This number of independent samples came as a result of a 
monochromatic transmission. Suppose the scatterometer 
transmission frequency is increased by Af = 5 MHz from 
pulse to pulse for twenty consecutive pulses and then re- 
peated. This corresponds to using M = 20 m Fig. IV- 1. 
The number of independent samples using this frequency 
agility would be twenty times the number arrived at 
using the monochromatic signal For the above mentioned 
transmission and reception specifications, at least 
20 x 18 = 360 independent samples would be received. 

This is dependent upon the requirement that the extreme 
conditions used to calculate Af c will be applied to the 
scatterome ter. 

In using a frequency variation, the equation for the 
backs catter coefficient will have an inherent error, as 
stated before, due to its calculation by the scatterome- 
ter’s data normalizer being dependent upon a specific 
wavelength or frequency value. The backs catter coeffi- 
cient, Ooj can be considered as a function of either 
^■ 3 or f 3 . The variance m o 0 is calculated as follows 
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t £tr 0 - I ~(i + "f ^) 3 

: I _ ( , + nr) 

= .01 

in db , this will be 

£02 = - I 0 log LOI 

= ±. 0 H db 

Thus, for the case cited, m causing an inherent db 
variance of ±.04 db m the calculated value of Oo , the 
number of independent samples are increased from 18 to 
20 x 18 or 360. The use of frequency agility must then 
be considered to be advantageous from the standpoint of 
the ability to increase the number of independent samples 
with only a small inherent error m o 0 . 
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CHAPTER V 
CONCLUSION 

A system analysis was performed to determine the 
effectiveness and accuracy of measuring backs catter co- 
efficient for a proposed scanning radar scatterometer. 

The scatterometer system was examined m detail by in- 
dividual subassemblies transmission, reception, signal 
reduction and logic, recording, and CRT display The 
system calibration and klystron operations were presented 
and analyzed m regards to implementation effectiveness 
and to transmission frequency drift respectively The 
possibility of improper CRT display of the obtained back- 
scatter coefficient values was presented, and a method 
for correcting the problem is suggested. 

Although the problems of transmission frequency 
drift and improper CRT display illumination were shown to 
be either mconsequencial or easily resolved, the problem 
of calibration presented some difficulty By mixing the 
calibration and intermediate frequency signals of equal 
amplitude with the return signal spectrum, improper power 
distribution was shown to occur. Equipment capable of 
handling a calibration signal of amplitude much smaller 
than the intermediate frequency signal amplitude could 
possibly solve this situation, but efforts must be made 



123 


to assure that the calibration signal amplitude is 
greater than any system noise level An alternate effec- 
tive method for injecting the calibration signal into the 
scatterometer system was presented. 

The accuracy of the mean determi nance of a back- 
scatter signal using a half-wave linear detector was 
analyzed. The analysis was conducted such that all re- 
sults are applicable to the proposed scatterometer sys- 
tem. Since the accuracy of a stochastic signal process 
must be defined m terms of signal independence, an auto- 
correlation function was derived which describes signal 
independence and is applicable to the scatterometer sys- 
tem The mean value of the input signal was assumed to 
be zero m order to derive the autocorrelation function. 
This being a necessary, but incorrect, assumption, strin- 
gent specifications as to what constitutes signal inde- 
pendence were defined for the autocorrelation function 
The scatterometer reception of the Rayleigh distributed 
voltage return signal was shown to contain 18 independent 
samples during the 0 2 second integration time for a 
doppler bandwidth of 280 Hz. Averaging 18 independent 
samples gave a level of confidence of 73% that the back- 
scatter coefficient is correct to within ±1 0 db of the 
true value These results are considered inadequate, an 
adequate level of confidence being considered greater 
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than 90%. It was shown that adequate results can be ob- 
tained by averaging at least 49 independent samples 
This increase in independent samples can be obtained by 
either increasing the doppler bandwidth or the viewing 
time for a terrain cell. Since the minimum doppler fre- 
quency used m defining a terrain cell size is utilized 
to determine the number of independent samples, the ac- 
curacy of the calculated backscatter coefficient for this 
cell can be considered the worst case. 

The advantages of implementing frequency agility m 
the transmission section of the s catterome ter were de- 
rived utilizing a specified frequency autocorrelation 
function. The number of independent sample was shown to 
greatly increase by utilizing frequency agility The use 
of frequency agility m the s catterome ter was deemed ad- 
vantageous from the standpoint of improving the accuracy 
of the backs catter coefficient de terminance . The main 
disadvantage incurred utilizing frequency agility is that 
of increased hardware complexity. 

In order to increase the s catterome ter accuracy 
without increasing the original integration time and ter- 
rain cell size, the effects of utilizing a frequency 
agility transmission system were investigated. Using a 
specified frequency autocorrelation, the implemented fre- 
quency change necessary to gam pulse- to-pulse signal 
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independence was found to be about 5 MHz An example, 
using the proposed scat terome ter specifications, showed 
that the utilization of frequency agility was very ad- 
vantageous from the standpoint of accuracy 

The frequency autocorrelation utilized is a function 
of target length and not target roughness A frequency 
autocorrelation function dependent upon the surface de- 
correlation distance has been attempted [33] Further 
work dealing with frequency autocorrelation functions for 
rough surfaces is justified and necessary before any lm- 

l 

plementation of frequency agility m the field of scan- 
ning radar scatterometers 

It is believed that future derivations concerning 
frequency autocorrelation functions should relate the 
critical frequency change and the surface decorrelation 
distance as directly proportional Thus, the rougher the 
surface, the smaller the frequency change Exactly how 
this function will appear is as yet uncertain, but its 
derivation could be a major factor m the design of fu- 
ture s catterometer systems 



126 


BIBLIOGRAPHY 


[1] Useful Applications of Earth - Oriented Satellites , 
National Academy of Sciences, Washington, D. C , 

1969. 

[2] R. K Moore, "Radar as a Sensor," University of 
Kansas, Lawrence, CRES Rept, 61-7, 1966 

[3] P Beckmann and A Spizzichmo, The Scattering of 
Electromagnetic Waves from Rough Surfaces New 
York The MacMillan Co., 1963, p. 119 

[4] J. Rouse, Jr., H. MacDonald, and W Waite, "Geo- 
science Applications of Radar Sensors," IEEE Trans - 
actions on Geoscience Electronics , vol. GE-7, 

pp 2-19, January 1969. 

[5] R. K Moore, "Radar Scatterometry - An Active Remote 
Sensing Tool," University of Kansas, Lawrence, CRES 
Rept 61-11, 1966 

[6] J W Rouse, Jr , "Arctic Ice Type Identification 
Radar," Proc, IEEE , vol. 57, pp 605-611, April 1969. 

[7] R W Newton, "Radar Scatterometer Data Analysis 
Sea State," Texas A§M University, College Station, 
Tech. Rept. RSC-09, 1970. 

[8] N Gurnard and J Daley, "An Experimental Study of 
a Sea Clutter Model," Proc IEEE, vol. 58, pp. 543- 
550, April 1970. 

[9] J W Rouse, Jr., "Radar Scatterometer Data Anal- 
ysis," Texas A§M University, College Station, 

Tech Rept RSC-01, 1969 

[10] J. W Rouse, Jr , "The Frequency Dependence of Back- 
scatter from Rough Surfaces," University of Kansas, 
Lawrence, CRES Rept. 133-4, 1968 

[11] W Ament, F MacDonald, and R Schewbridge, "Radar 
Terrain Reflections for Several Polarizations and 
Frequencies," Naval Research Laboratory, Unpublished 
Rept , 1959. 



127 


[12] R Cosgriff, W. Peake, and R Talor, "Terrain Scat- 
tering Properties for Sensor System Design," Ohio 
State University Engineering Experiment Station, 
Columbus, Bulletin No 29(3), 1960 

[13] J Lundien, "Terrain Analysis by Electromagnetic 
Means," U S. Army Engineer Waterways Experiment 
Station, Tech. Rept No. 3-693, 1966 

[14] D Earing, "Target Signature Analysis Center Data 
Compilation, Second Supplement," University of 
Michigan, Institute of Science and Technology, Ann 
Arbor, 1967. 

[15] W Birkemeier and N. Wallace, "Radar Tracking Ac- 
curacy Improvement by Means of Pulse- to -Pulse Fre- 
quency Modulation," Trans AIEEE, vol 81, pp 571- 
575, 1962 

[16] R E. Chapman, "Systems Specifications - Model 719 
Scanning Scat terome ter , " Ryan Aeronautical Co , 

San Diego, Internal Report, 1969 

[17] T K Ishu, Microwave Engineering . New York 
Ronald Press Co., 1966, p 231. 

[18] R. E Chapman, op cit 


[19] G W Dalke, J. E. Rathke , J W Rouse, Jr., 
"Initial Design Considerations for the Multi-Color 
Display," University of Kansas, Lawrence, CRES 
Tech. Memo 61-5, 1964. 

[20] A Papoulis, Probability, Random Variables, and 
Stochastic Processes. New York McGraw Hill, 
1965, p 100-104. 

[21] R. E Chapman, op . cit . 

[22] W Davenport, Jr and W Root, Random Signals and 
Noise . New York McGraw Hill, 1958, p 271 

[23] Ibid. , p 104 

[24] A. Papoulis, op cit , p. 266 

[25] A Papoulis, op. cit. , pp 266-268. 



128 


[26] 

[27] 

[ 28 ] 

[29] 

[30] 

[31] 

[32] 

[33] 


R. K. Moore, "Radar Scatterometry - An Active Remote 
Sensing Tool," ibid. 


A. 

Papoulis , 

op 

Clt. , 

p. 260 

A. 

Papoulis , 


Clt. , 

p 94. 

A. 

Papoulis , 

°£^ 

Clt. , 

p . 64 . 

A. 

Papoulis , 

op. 

Clt. , 

p. 129 


I. Katz, "Wavelength Dependence of the Radar Re- 
flectivity of the Earth and Moon," J. Geophysical 
Research, vol. 71, pp. 361-366, Jan. 15, 1966 

W. Birkemier and M. Wallace, op, cit. 


R. Moore, W. Waite, J. Rouse, Jr., "Panchromatic 
and Polyp an chromatic Radar," Proc. IEEE, vol. 57, 
pp. 590-593, April 1969. 



129 


VITA 

Richard Harrison Arnold 

|2BI the son of Mr. and Mrs Howard 

Harrison Arnold He received his elementary education at 
several schools and graduated from Lafayette High School, 
Lafayette, Louisiana, m May, 1965. 

Mr. Arnold then attended the University of South- 
western Louisiana at Lafayette. He graduated with honor 
m May, 1969, having become a member of Eta Kappa Nu, 

Tau Beta Pi, and IEEE, and received the Bachelor of 
Science Degree m Electrical Engineering. During the 
1968 summer secession, he was employed as a Student 
Change Engineer at The Boeing Co m New Orleans working 
with the Apollo Saturn V Boosters 

In June, 1969, Mr Arnold entered Texas A§M Univer- 
sity and completed the requirements for a Master of 
Science Degree m Electrical Engineering m August, 1970 
During his tenure at Texas A§M University, he was em- 
ployed as a Research Assistant by the Remote Sensing 
Center 

Mr Arnold is married to the former Clara Elizabeth 
Rountree and has a daughter Grace Elizabeth He can be 
reached at P 0. Box 468, Pmeville, Louisiana, 71360 

The typist for this thesis was Mrs. James Ronald Bush 





The REMOTE SENSING CENTER was established by authority of the Board of Directors of 
the Texas A&M University System on February 27 1968 The CENTER is a consortium of four 
colleges of the University, Agriculture, Engineering Geosciences and Science This unique 
organization concentrates on the development and utilization of remote sensing techniques and 
technology for a broad range of applications to the betterment of mankind. 





